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Method of Identification of Non-primary Events in Seismic Data 

This invention relates to seismic signal processing generally and specifically 
relates to the identification and attenuation of non-primary events in seismic data. 

hi seismic data processing, identification and attenuation of non-primary events 
(e.g. multiples, ghosts, etc.) is desirable. One specific example of such efforts is the so- 
called multiple diffraction attenuation (MDA), which is a processing technique that has 
been designed by PGS Data Processing to attenuate high-energy reverberations of near- 
surface diffractions. Such reverberations are commonly seen on deep-water datasets with 
complex, shallow sub-seabed topography such as those recently acquired by PGS 
offshore Nigeria and in the South Rockall Trough. 

Conventional MDA processing is based on the assumption that, at greater depths, 
multiples tend to have significantly higher amplitudes at higher frequencies than 
primaries do. The reason for this is greater primary absorption that is caused by (1) 
increased wave-field spreading (longer travel-paths) and (2) higher absorption. Figure 1 
illustrates the conventional procedure. In the high frequency constituent of the data, high 
amplitudes HA are identified by comparing amplitudes to the rms amplitude in a user- 
defined window. The boundaries of the window are dependant upon the specific shapes 
of the events expected bases on a priori knowledge and are within the skill of those in the 
art. For example, when amplitudes are higher than, say, two times the rms amplitude, 
then these amplitudes are zeroed. After stacking the high and low frequency constituents, 
the conventional MDA result is obtained. Although this method has proven to be an 
effective noise removal tool, the major disadvantage of the method is that it cannot 
separate high-amplitude primaries from high-amplitude multiples. As a result, high- 
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amplitude primaries may be partly muted, as indicated by the discontinuous primaries P. 
Clients often recognize this disadvantage; and, as a result, the method is often not 
applied. 

Therefore, there is a continuing need for a method of identification non-primary 
events, and it is an object of the present invention to provide for that need. 
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SUMMARY OF THE INVENTION 

According to one aspect of the present invention, a method is provided for 
identifying primary seismic events in seismic data, the method comprising: 

applying a coherency filter to the seismic data; 

sorting the data according to an event characteristic having a tendency to separate 
primary from non-primary events; and 

selectively attenuating according to the characteristic and the coherency of the 

events. 

According to a more specific example of the present aspect, the coherency filter is 
applied in windows, the characteristic comprises amplitude in a limited range of 
frequencies, and the attenuation comprises reduction of amplitude (e.g. muting). 

According to another broad aspect of the invention a system is provides for 
identifying primary seismic events in seismic data, the method comprising: 

means for applying a coherency filter to the seismic data; 

means for sorting the data according to an event characteristic having a tendency 
to separate primary from non-primary events; and 

means for selectively attenuating according to the characteristic and the coherency 
of the events. 

According to still another aspect of the present invention, there is provided a 
method of identification of primary events in seismic data, the method comprising: 

sorting the data by frequency wherein at least some non-primary events are 
separated from primary events, wherein a frequency-sorted gather of data results; 
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attenuating in the frequency-sorted gather amplitudes above a pre-selected base 
amplitude, wherein attenuated amplitudes result; 

applying a coherency filter to the events, wherein coherent events are identified; 

replacing with amplitudes from the coherent events attenuated amplitudes in the 
frequency-sorted gather corresponding to the coherent events. 

According to a more specific example of the present aspect, the attenuating 
comprises reducing amplitude. According to an alternative example, the attenuating 
comprises muting. 

According to another aspect of the present invention, there is provided a method 
of identification of primary events in seismic data, the method comprising: 

sorting the data by frequency wherein at least some non-primary events are 
separated from primary events, wherein a frequency-sorted gather of data results; 

applying a coherency filter to the events, wherein coherent events are identified; 

attenuating in the frequency-sorted gather amplitudes above a pre-selected base 
amplitude which are not associated with the coherent events, wherein attenuated 
amplitudes result. 

According to a more specific example of the present aspect, the attenuating 
comprises reducing amplitude. According to an alternative example, the attenuating 
comprises muting. 

According to yet another aspect of the present invention, there is provided a 
system of identification of primary events in seismic data, the method comprising: 

means for sorting the data by frequency wherein at least some non-primary events 
are separated from primary events, wherein a frequency-sorted gather of data results; 
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means for attenuating in the frequency-sorted gather amplitudes above a pre- 
selected base amplitude, wherein attenuated amplitudes result; 

means for applying a coherency filter to the events, wherein coherent events are 
identified; 

means for replacing with amplitudes from the coherent events attenuated 
amplitudes in the frequency-sorted gather corresponding to the coherent events. 

According to a more specific example of the present aspect, the means for 
attenuating comprises means for reducing amplitudes. According to another example, the 
means for attenuating comprises means for muting. 

According to still a further aspect of the present invention, a system of 
identification of primary events in seismic data is provided, the system comprising: 

means for sorting the data by frequency wherein at least some non-primary events 
are separated from primary events, wherein a frequency-sorted gather of data results; 

means for applying a coherency filter to the events, wherein coherent events are 
identified; 

means for attenuating in the frequency-sorted gather amplitudes above a pre- 
selected base amplitude which are not associated with the coherent events, wherein 
attenuated amplitudes result. 

According to a more specific example of the present aspect, the means for 
attenuating comprises means for reducing amplitude. According to an alternative 
example, the means for attenuating comprises means for muting. 
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BRIEF DESCRIPTION OF THE DRAWINGS: 

Figure 1 is a plot of a set of seismic data used in accordance with and example 
embodiment of the present invention. 

Figure 2a is a plot of a set of seismic data used in accordance with and example 
embodiment of the present invention. 

Figure 2a is a plot of a set of seismic data used in accordance with and example 
embodiment of the present invention. 

Figure 3a - 3e is a plot of gathers of seismic data used to illustrate examples of 

the present invention. 

Figure 4 is a block diagram of an example embodiment of the present invention. 
Figure 5 is a block diagram of an example embodiment of the present invention. 
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DESCRIPTION OF EXAMPLE EMBODIMENTS OF THE INVENTION: 

Referring now to Figure 2a, an example of the present invention will be illustrated 
in which coherency filters are applied in windows W surrounding high amplitude events 
HA in the high frequency content of seismic survey data. Various forms of coherency 
filters are useful according to the invention, and their design and selection will be within 
the skill of those in the art upon review of the present description. In the illustrated 
example, NMO-corrected CDP gathers are seen, although other gathers that will occur to 
those of skill in the art will also suffice. After application of the coherency filter, a 
smoothed CDP gather is obtained in which non-coherent events are significantly 
attenuated. Primary events are more coherent that non-primaries. 

The seismic data are also subjected to processing in which amplitudes HA over a 
pre-determined value in the higher frequency components of the data are attenuated. The 
amplitudes HA for attenuation are determined, as will be understood by those of skill in 
the art, by examination of the frequency content and amplitudes of specific data sets. In 
some embodiments, the amplitudes HA for attenuation comprise rms amplitudes, while in 
other embodiments, the amplitudes HA for attenuation comprise peak amplitudes. Still 
other amplitude measurements within the scope of the present invention will occur to 
those of skill in the art. 

Likewise, the distinction between those frequencies in which significant 
amplitude differences exist between primaries and non-primary events is survey specific 
and is selected based on examination of the particular data. In an alternative 
embodiment, however, the pre-determined value of the amplitude is not based on data 
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examination, but rather it is chosen by the experience of the processor or pre-set in the 
processing software on which the method of the present invention is run. 

Referring now to Figure 2b, the attenuated samples high-frequency amplitude 
attenuation are replaced by the corresponding samples from the coherency filter method. 
As a result, high-amplitude, coherent primaries P are maintained, whereas non-coherent 
multiples are strongly attenuated. 

Referring now to Figures 3a-e, the application of an example method of the 
present invention is seen. The first gather (Fig. 3a) shows the high-frequency constituent 
of an arbitrary CDP gather after NMO correction. The second gather (Fig. 3b) shows the 
same gather, after application of coherency filtering. The third gather (Fig. 3c) shows the 
same gather after application of conventional MDA. Note the strong attenuation of the 
high-amplitude primaries. The fourth gather (Fig. 3d) shows the result after application 
of the present invention. Note the significant reduction of incoherent high-amplitude 
events, as well as the preservation of the coherent events that are believed to be primaries. 
Finally, the fifth gather (Fig. 3e) shows the difference plot between the present invention 
and conventional MDA. 

The above method is carried out according to various systems. In one such 
system, illustrated in Figure 4, means 40 for sorting the data D by frequency separates at 
least some non-primary events (NPE) from primary events (PE), wherein a frequency- 
sorted gather of data (FSG) results. The means 40 is in communication with means 44 
for attenuating in the frequency-sorted gather (FSG) amplitudes above a pre-selected base 
amplitude, wherein a new frequency-sorted gather with selected attenuated amplitudes 
(FSGAA) results. In turn, means 40 is in communication with means 46 for applying a 
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coherency filter (CF) to the events, wherein coherent events (CE) are identified. Finally, 
means 44 and means 46 are connected to means 48 for replacing with amplitudes from 
the coherent events attenuated amplitudes in the frequency-sorted gather corresponding to 
the coherent events. The result is a frequency-sorted gather with attenuated non-primary 
events with substantially unattenuated primary events. 

According to still a further aspect of the present invention, an alternative system 
of identification of primary events in seismic data is provided, the system comprising 
means 40 for sorting the data by frequency wherein at least some non-primary events 
(NPE) are separated from primary events (PE), wherein a frequency-sorted gather (FSG) 
of data results. The means 40 is connected to means 46 for applying a coherency filter 
(CF) to the events, wherein coherent events (CE) are identified. In this embodiment, 
means 44 attenuates in the frequency-sorted gather amplitudes above a pre-selected base 
amplitude which are not associated with the coherent events, wherein a frequency-sorted 
gather with attenuated non-primary events but substantially unattenuated primary events 
results. 

In one more specific example of the systems of Figures 4 and 5, the means 44 for 
attenuating comprises means for reducing amplitudes. According to another example, the 
means for attenuating 44 comprises means for muting. 

As will be understood by those of skill in the art, the various means 40 - 46 
comprise, in some embodiments, hard-wired circuitry (either analog, digital, or some 
combination of the two) or, alternatively, computer hardware, programmed as will occur 
to those of skill in the art to perform the functions described. In still other embodiments, 
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the means 40 - 48 comprise objects and classes of code which are connected in a 
computer as will occur to those of skill in the art. 

The above embodiments are given by way of example, only. Other embodiments 
will occur to those of skill in the art that do not depart from the spirit of the present 
invention. 
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